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Precursor of Color Superconductivity in Hot Quark Matter 
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We investigate possible precursory phenomena of color superconductivity in quark matter at finite 
temperature T with use of a low-energy effective theory of QCD. It is found that the fluctuating pair 
field exists with a prominent strength even well above the critical temperature Tc- We show that the 
collective pair field has a complex energy located in the second Riemann sheet, which approaches 
the origin as T is lowered to Tc. We discuss the possible relevance of the precursor to the observables 
to be detected in heavy ion collisions. 
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It is an intriguing subject to explore how the QCD 
vacuum changes in hot and dense hadronic matter j^, ^. 
The recent lattice QCD simulations |^ show that the 
QCD vacuum undergoes a phase transition to a chirally 
restored and deconfincd phase at as low temperature (T) 
as T ^ 170 MeV. Although it is still a great challenge to 
perform a Monte Carlo simulation with a finite chemical 
potential /i on the lattice QCD with SU(3) color group, 
we believe that a chirally restored and deconfincd system 
is realized in dense hadronic matter with large fj,. 

It is to be noticed that a system with a finite fi at 
low T can have a Fermi surface and the system may be 
described as a Fermi liquid composed of quarks. In ac- 
cordance with the validity of the diquark-quark picture 
of baryons[Q, the quark-quark interaction is attractive in 
some specific channels. Then the existence of the Fermi 
surface gives rise to a Cooper instability with respect to 
the formation of the diquark or Cooper pair in the most 
attractive channel, and the system is rearranged to a su- 
perconducting phase where the color-gauge symmetry is 
dynamically broken^. Such a color superconductivity 
(CSC) has recently acquired a renewed interest, because 
the resulting gap A was shown to be as large as 100 
MeV[^, 0. Furthermore the color and flavor degrees of 
freedom of quarks give a fantastically rich structure in 
dense quark matter with CSCj|]: Such a quark matter 
may be realized in the deep interior of neutron stars (NS) , 
and characteristic phenomena observed for NS such as 
the glitch phenomena might be attributed to the exis- 
tence of the CSC in the NS. 

Can CSC be relevant to experiments in the labora- 
tory on the earth? The experiment using the available 
facilities, such as RHIC using high energy heavy ion colli- 
sions seems unfortunately not fit for the problem, because 
the matter produced by RHIC is almost baryon-free with 
much higher T than the critical temperature Tc ^ 50-60 
MeV of CSC. The purpose of the present paper is to ex- 



plore possible precursory phenomena of CSC at T > Tc, 
i.e., in the Wigner phase and show that the fiuctuation 
of the pair field is significant at T well above T^, which 
thereby might affect observables to be detected by the 
heavy ion collisions with large baryon stopping 

Since we are interested in relatively low T and p region, 
we adopt a low-energy effective theory of QCD as was 
done in |l^ where an instanton-induced interaction is 
employed. We remark here that a simplified version of 
the instanton-induced interaction in the two-flavor case 
takes the form of Nambu-Jona-Lasinio (NJL) modeljlj] 
which can be also derived by a Fierz transformation of 
the one-gluon exchange interaction with heavy-gluon ap- 
proximation; see |jl|, ^, |l3|, Indeed it was shown later 
|lT| that the physical content given in ||l^ can be nicely 
reproduced by the simple NJL model. Therefore we sim- 
ply take the following NJL model as an effective theory. 
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where V''^ = Cip'^ , with C — Z7^7° being the charge con- 
jugation operator. Here we have confined ourselves in 
the two-fiavor case in the chiral limit; T2 and A2 are the 
second component of the Pauli and Gell-Mann matrices 
representing the flavor SU(2) and color SU(3), respec- 
tively. The scalar coupling constant G5 and the three 
dimensional momentum cutoff A are determined so that 
the pion decay constant /^r = 93 MeV and the chiral 
condensate {-ipip) — (— 250MeV)'^ are reproduced in the 
chiral limit[|l3|; Gs = S.OlCeV"^ and A = 650 MeV. As 
for the diquark coupling constant Gc , we shall fix Gc at 
3.11GeV"^ in this paper. We have checked that the 
qualitative features of the results do not change for Gc 
varied in the range used in the literature j|, [TT| . 

The phase diagram determined in the mean-field ap- 
proximation with is shown in Fig.l: The dashed line 
denotes the critical line for a second order phase transi- 
tion and the solid line first order phase transition. The 
critical temperature of the chiral transition at /i = is 
T = 185MeV, while the critical chemical potential of 
chiral-CSC transition at T = is /i = 316MeV. The 
tricritical point where the order of the chiral transition 
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FIG. 1: The calculated phase diagram in T-fi plane in our 
model. The Solid and dashed line denote the critical line of a 
first and second order phase transition, respectively. 



changes from second to first order is located at about 
(T,/z) = (70, 280)MeV. This phase diagram essentially 
coincides with those given in |l^ . 

Given the phase diagram for the equilibrium state, we 
now proceed to examine fluctuations of the color Cooper 
pairs in the normal phase a,t T > Tc. We shall obtain 
possible collective states and their strength functions in 
the hnear response theory. We suppose that the follow- 
ing external pair field is applied to the system, A*^(a:) = 
^exS**' ''fi"'""'^, which is tantamount to adding the exter- 
nal Hamiltonian H^x = / rfxA*^(a;)('!/''^«75T2A2'(/') + h.c. 
Then a color Cooper-pair field is dynamically induced, 
as given by A*„^(x, r) = -2Gc(^(x, T)T2A2n5V''^(x, r)), 
where the bracket denotes a statistical average. Using the 



anomalous Green function defined by -F'l'(x, r; x', r') = 
-(T.r{V''^(x,r)V;(x',r')}), the induced pair field is ex- 
pressed as 

A*„,(x,t) = -2GcTr T2A2Z75Ft(x, r; x, r+), (2) 

where Tr denotes the trace taken over the color, flavor 
and spinor indices. Ft(x, r; x', r') satisfies the following 
differential equation 

[-7°|:-i^o-A*7°]i^t(x,T;x',T') 

= ^tot(x,'r)T2A2i75G(x,T;x',T'), (3) 

where -ftTo = -if- V + A/ and M = -2Gs(-0(x, r)?A(x, t)) 
is the dynamical mass. Here, A^^^{x,t) — A*j.(x, t) -|- 
A*j^(x, r) is the total pair field to induce the pair 
field: Notice the self-consistent nature of the prob- 
lem. G(x, r;x',r') is the usual non-anomalous Matsub- 
ara Green function. To see a linear response of the 
system, one may replace G(x, r;x',r') in Eq.(0) with 
a free Matsubara Green function G'-°-'(x, t; x', r^; thus 
we have Ft(x, r; x', r+) = (-7"!: - Ko - ^7°)"' x 
A(Q((x, r)r2A2Z75G^°^(x,r;x, T+). Inserting this expres- 
sion into Eq.(H) and performing the Fourier transforma- 
tion, we have 

A*„rf(k,w„) = -GcA*„,(k,c^„)Q(k,^„), (4) 

where 
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(5) 



Here, E^ = ^/WTM^ and J^{E) = [exp(/3(£; ± ii)) + 
with (3 = T~^. The first term in the brace on the 
r.h.s. of Eq.(^ is originated from the process A^^j-l-g q 
and its inverse process q — *■ l^tot + 9- The second and 
third terms can be interpreted similarly in terms of ki- 
netic processes. Inserting A(Q((k, aj„) = A*^(k, cj„) + 
A*^^(k, w„) into Eq.(U), we obtain the induced pair field 
in the linear response to the external field in the imagi- 
nary time as 

^-d(k,^„) = T^^^^fe^AL ^ I?(k,za;„)AL. (6) 



This is a central result in the present work, on which 
the rest of discussions is based. The proportionality 
constant in Eq.(^ converted in the real-time domain is 
called the response function or retarded Green's function 
£'^(k, cj), which is obtained by the analytic continuation 
of I?(k, bj) to the upper half of the complex energy plane 
as D^{k, Lo) = ^(k, bJ + ir]). 

Setting M — 0, let us examine precursory collective 
modes above the CSC phase when approaching the crit- 
ical point as shown by the arrow in Fig. 1. Then we 
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notice the following equality holds 

1 + Gc 2(0, 0)1^,^3,^ =0, 



(7) 



owing to the self-consistency condition for the diquark 
condensate at T — T^. This equality is the origin to lead 
to the various precritical phenomena for CSC , as we will 
see shortly. Here it should be stressed that although the 
following discussion is seemingly relied heavily on the fact 
that phase transition is second order, the nature of the 
collective mode to be discussed below will not be altered 
even when the phase transition is weak first order. 

The spectral function which is the measure of the ex- 
citation strength of the fluctuating pair field is given by 



p(k,w) = --ImD^(k,w). 



(8) 



The temperature dependence of the spectral function at 
/i = 400 MeV and zero momentum transfer (k = 0) is 
given in Fig. 2 (a). One can see a prominent peak moving 
toward the origin grows as the temperature is lowered 
toward Tc as indicated by the arrow in Fig.l. In the 
experimental point of view, it is interesting that the peak 
survives at T even well above Tc with e = (T — Tc)/Tc ~ 
0.2. This means that the precritical region of CSC is 
one to two order larger in the unit of Tc than that in 
electric superconductors. Of course, this is related to the 
difference of the ratio ujd/Ef and A/Ep where lud is 
the Debye frequency and Ep is the Fermi energy; NJL 
model is similar to the attractive Hubbard model than to 
BCS modclllll. The result for fi = 500 MeV is shown in 
Fig. 2(b); one sees that although the growth of the peak 
becomes relatively moderate, the qualitative features do 
not change. 
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FIG. 2: The spectral function for the pair field at T > Tc 
with e= [T - Tc)/Tc = 0.05,0.1,0.2 and 0.5 at ^ = 400MeV 
(a) and fj. = 500MeV (b). 



The existence of the peak with the narrow width sug- 
gests an existence of a well-defined collective mode as an 
elementary excitation in the Wigner phase. The exis- 
tence and the dispersion relation uj = Lu{k) of a collective 
excitation in a isotropic system is examined by searching 
possible poles of the response function, i.e., solving the 
equation 1 -I- GcQ{k,uj) = 0. Here Q{k,uj) is the an- 
alytic continuation of Q{k,uj„) to the upper half plane. 



The existence of a pole ui = uj{k) for given k means that 
the pair field Ai„(i(fc, a;(fc)) can be created even with an 
infinitesimal external field A*^ as seen from Eq. , which 
physically means that the system admits spontaneous ex- 
citation of a collective mode with the dispersion relation 
UJ = uj{k). Since the energy of the collective mode is 
found to be the lowest at A: = 0, let us obtain uj{0) = z. 
One should notice that the gap equation Eq.(0) implies 
that z = is a solution at T = Tc, i.e., there exists a zero 
mode at the critical point. The singular behavior seen 
in the spectral function near the critical point in Fig. 2 is 
caused by the presence of the zero mode. Conversely, one 
may take the emergence of the zero mode as the critical 
condition for the color superconductivity. This condition 
is actually long known in the condensed matter physics 
and called Thouless criterion JTsl. 

Notice that the pole should be located in the lower half 
plane otherwise the system is unstable for the creation of 
the collective mode. Because of a cut along the real z 
axis, the analytic continuation of (5(0, z) to the lower 
half plane leads to the second Riemann sheet; 



2-2^z(^ -I- -)2 tanh — ] , 



(9) 



where P denotes the principal value and the term in the 
second line is the additional term for z to be in the second 
sheet. 

A numerical calculation shows that there indeed exists 
a pole in the lower half plane for T > Tc- Fig. 3 shows how 
the pole moves as T is lowered toward Tc for = 400 and 
500 MeV. One sees that the pole approaches the origin as 
T is lowered toward Tc- Such a mode whose energy tends 
to vanish as the system approaches the critical point of 
a phase transition is called a soft mode0|. Thus we have 
found the pole for the soft mode for CSC in hot quark 
matter. As far as we know, this is the first time when the 
identification has been done of the softening pole in the 
complex energy plane for the superconductivity including 
the condensed matter physics p5[| . The significance of the 
existence of such a soft mode at fc = lies in the fact that 
it causes a long-range correlations for the color Cooper 
pairs in the Wigner phase, which will leads to a singu- 
lar behavior of various observables such as the transport 
coefficients, as will be discussed later. We remark that 
such a soft mode in hot and dense quark matter was first 
discussed in [ pTf in relation to the chiral transition. 

Another characteristic feature of the soft mode is that 
the absolute value of the imaginary part uji of the pole 
is larger than that of the real part lOj. for both /i. This 
feature may have an important implication for an effec- 
tive equation describing the dynamical phase transition 
for CSC: The dynamical behavior of the order parameter 
near Tc is well described by a non-linear diffusion equa- 
tion which is known as the time-dependent Ginzburg- 
Landau (TDGL) equation (l8|. We however remark here 
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that by setting k = 0, the so called Landau damping 
due to the particle-hole excitation is not incorporated in 
the present calculation. Fig. 3 shows that the ratio coi/ujr 
becomes larger for the larger fj,. This is plausible be- 
cause the larger the Fermi energy, the larger the density 
of states near the Fermi surface, which implies that the 
volume of the phase space for the decay of the pair field 
increases with fi at low T. 
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FIG. 3: The pole of the precursory soft mode for /i = 400 
MeV and 500MeV. The numbers attached to dots denote e — 
(T - TJ/T,. 



It is interesting that the so called pseudogap in the 
density of states observed in high-T^ superconductors 
(HTSC) as a characteristic precursory phenomenon! 19 
might be attributed to fluctuations of the Cooper 
pairs Thus the present results strongly suggests such 
a pseudogap phase exists also for CSC in a rather broad 
region of T above Tc. To establish the pseudogap phase, 
one must calculate the density of states from the single 
particle Green's function Jl6|, |20|] , which we hope to report 
in a future publication. Experimentally, a measurement 
of the single-particle quark propagator is useful for ob- 
taining the information on the density of states for the 
single quark states, as was done for HTSC by a kind of 



the quasi-elastic scattering by photon, called the angle- 
resolved photo-emission spectroscopy (ARPES)|l9). 

The fluctuation of the pair field in metals above Tc 
affects several transport coefficients like the electric con- 
ductivity (EC), the phonon absorption coefficient and the 
specific heat|gl[. An anomalous excess of the electric 
conductivity Tc is known as paraconductivity (PC)pll. 
The direct analogy thus leads to an expectation of the 
enhancement of color conductivity above Tc in our case, 
which we call color-paraconductivity (CPC): In fact, the 
PC for the metal superconductivity comes from the fact 
that the effective dynamics is given by TDGL equation 
giving rise to a zero mode at momentum k — at 
T = Tc ||2l|, which is just the behavior extracted from 
the behavior of our soft mode as shown in Fig. 3. 

More interesting is that the electric current can also 
have an anomalous enhancement in the Wigner phase 
due to the color pair- field fluctuations: This is because 
the propagator of photons (and also gluons) at T > Tc 
is modified by the fluctuating pair-field as in the usual 
superconductors [pi] . This implies that photons and 
hence dilcptons from the system can carry some infor- 
mation of the fluctuation of the color pair field. 

In short, we have explored precursory phenomena for 
the color superconductivity (CSC) in quark matter at 
T =/= 0; the fluctuating pair field exists with a prominent 
strength even well above Tc. We have also discussed the 
possible observability of CSC through the fluctuations 
by heavy ion collisions with large baryon stopping using 
an analogy with the precritical phenomena known in the 
metal and the high- Tc superconductors. 
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